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ABSTRACT 


The reliability of the structure can be calculated by combining the 
non-uniform stress distribution of a structure from finite element analysis 
with the statistical strength theory under two dimensional non- 
homogeneous, uniaxial stress limitations. The specimen models with and 
without notch are the sample structures to illustrate the calculations. The 


Statistical strength of the structure is cast in the standard Weibull form 


characterized by the structural scale parameter, Br and shape parameter, 


ox which are functions of the material scale and shape parameters («x and 


8.) Therefore, this thesis demonstrated that the scale parameter of the 
Structure remains constant for different load magnitude, only the 


structural geometry and loading condition cause the change of material 
scale parameter, 8 to the structural scale parameter, Be. The results also 
shows that the degree of uniformity of stress within structural element 


and the linearity of the load magnitude and stress distribution within the 


Structure affect the accuracy for the calculation. 
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. INTRODUCTION 


A. RELIABILITY REQUIREMENT FOR COMPOSITE STRUCTURES 

Reliability theory became an independent scientific discipline at the 
beginning of the 1950’s under the influence of the rapid development of 
radio electronics, computer, devices, and rocketry. Modern radio electronic 
apparatus and digital mechanics consist of a very large quantity of 
components. If the failure of one component leads to the failure of the 
unit as whole, then evidently the possibility of high reliability (freedom 
from failures) will diminish rapidly as the unit becomes more complex. In 
this connection, problems arise on the prediction of the reliability of 
planned units, on the development of measures to increase the reliability, 
on a fundamental for reliability testing methods, etc. All these questions 
are the subject of reliability theory. 

The mechanical strength and stiffness of a structure are one of the 
aspects of reliability. The engineering design concept is not only solving 
the stresses and strains which originate in structures under various 
external imposed conditions but also trying to make them to be operated 
sufficiently reliably throughout their established utilization — time. 
Therefore, in the concluding stages of an engineering design structural 
mechanics inevitably comes into resolution with reliability theory. For 
critical application of composite as an aircraft structure, the study of 


reliability in composite structures is required. 


B. HISTORICAL BACKGROUND 

Bullock [ Ref. 1 ] used Weibull theory to predict the strength of the 
composite materials between different laboratory sample configurations. 
He investigated the strength ratios of composite materials in flexure and in 
tension. He assumed the shape parameter (the exponent characterize 
flaw-density that determines the scatter of strength of the material) was 
not changed when the geometry or size of the material was changed. He 
used the same shape parameter for the tensile test and flexure test which 
based on homogeneous state of stress and heterogeneous state of stress 
respectively. He found, however, based on the experimental results, the 
hypothesis that the shape parameter was constant was not realistic. 
Changes in shape parameter were in fact observed. Whitney and Knight, 1 
[ Ref. 2 ] used a statistical strength theory based on Weibull distribution to 
explain the difference between unidirectional tensile data generated from a 
flexure test and a standard tensile coupon. The result was shown a 
significantly larger variation in tensile strength versus flexural strength. 

Rosen [ Ref. 3 | presented a theoretical and experimental treatment of 
the failure of composites, consisting of a matrix stiffened by uniaxially 
oriented fibers, when subjected to a uniaxial tensile load parallel to the 
fiber direction. He observed that a portion of the fiber at each end was not 
fully effective in resisting the applied load because of the axial load was 
transmitted by shear through the matrix to adjacent fibers. Basically, the 
fibers failed as a result of a statistically distributed flaws or 
imperfections and composites failed as a result of a Statistical 


accumulation of such flaws over a givenregion. Therefore, the ineffective 


1S 


length of the fiber was introduced in order to consider the composite to be 
composed of a series of layers of this length. Then the segment of a fiber 
within a layer might be considered as a link in the chain that contribute the 
fiber. Each layer was then a bundle of such links, and the composite as 
whole can be modeled by series of such bundles. Weibull theory was used 
to define a statistical distribution of flaw or imperfections that result in 
fiber failure under applied stress. The statistical accumulation of such 
flaws within a composite material was demonstrated to be cause of 
composites failure. He concluded that the stress concentrations occurred 
In fiber adjacent to the fibers break were considered only to balance the 
effect of the variation of ineffective length which affect to the composite 
Strength. Knight, C.E. [ Ref. 4] used a finite element stress analysis and 
Statistical strength theory to assess the influence of the stress 
concentration on the ultimate strength calculation for composites from 
the Split-D Test. He observed that the geometry changed causes the 
nonlinearity and the boundary condition would be changed. 

The Probabilistic Statistic Failure of Composite Materials [ Ref. 5 | 
demonstrated how to develop a theory that combined the statistics of 
composite material failure with the orthotropic nature of composite 
materials. The effects of loading history and the probabilistic location of 
the failure for a composite can be accounted using this theory. In order to 
improve the efficiency of Composite designs for structural analysis, a 
theory that Combining the statistics of the experimental failure data with 
an orthotropic material failure response was developed. The analytical 


theory developed extent deterministic failure laws to include the 
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statistical scatter observed experimentally. Tsai and Wu [ Ref. 6 | are one 
of the investigators who described a general theory for the failure of 
composite materials as the strength tensor theory but determining is 
complicated by wide range in failure properties of composite materials and 
scatter in experimental results. The Weibull weakest link theory was 
developed and described a description of the spatial distribution of the 
failure locations with a sample calculation performed for the case of a 
uniformly loaded tensile specimen. Three point bending specimens were 
used to predict failure location by applying the basic weakest link theory. 
Four point bending specimens were used to illustrate the effect of loading 
history which were evaluated by applying numerical integration to the 
probability of failure density distribution curves. The compressive failure 
load was demonstrated that it is not deterministic failure and also having 
some scatter which is independent of tensile scatter so that the Weibull 
failure distribution for tension and compression are different. 

In ordinary structures having many stress components active, a 
decision on the probability for survival must be based on statistical 
descriptions that admit multiaxial stress. The development of combining 


weakest link theory with the strength tensor representation is required. 


C. OBJECTIVE AND SCOPE 

The driving motivation for this research is that many critical 
components of modern aircraft structure are made of composite materials. 
However, the reliability design is necessary for engineering designs 


requires not only the stress analysis but also statistical theory of failure 
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strength of composites, The latter motivated the investigation of 
statistical theory for composites failure strength where are summarized in 
the last section. Any aircraft structural components inevitability have the 
non-uniform stress in the neighborhood of edge-notches and holes. Size 
effect is the one that is importance for the composite study because of the 
difficulty and the economy that one could not experiment the several large 
samples to predict the reliability of the components. Even the 
experimental test samples can be made, the difficulty to identify the 
failure location is also impossible. Finite element method and the 
Statistical theory can be combined to explore the possibility of prediction 
both the reliability of the structure and the failure location. Basically, 
finite element method can be used to calculate the stress distribution that 
occurs in the structure under an arbitrary loading boundary condition. The 
Structure is divided into elements be sufficiently small that the stress is 
uniform. This solution of physically continuous system can then be 
replaced by the solution of the discrete system of the elements. The 
Statistic Weibull theory can be used to model composite materials 
composing of a series of bundles which are in turn composited series of 
fiber links. Another assumption for the statistical Weibull theory is the 
local stress in each link must be uniform. For the connection to the two 
dimensional finite element method, instead of using link, the element will 
be introduced. For one dimensional case, the link can be used but in two 
dimensional case, the element has a width. The operation of size effect 
formulation based on Weibull theory and finite element method will be 


presented in next chapter. 


The scope of this research limits to the treatment of statistical 
strength of composite structures which are two dimensional, subjected to 
nonhomogeneous uniaxial stress including size effect. The specimen 
models are plate with and without notch subjected to the tensile 
displacement loading case. The finite element code name NIKEZ2D [ Ref. 7 ] 
is used. NIKE2D is a two dimensional finite element program which is 
developed by Hallquist, J.O. [ Ref. 7 |. It can be used to calculate linear and 
non-linear finite element problems. In this thesis, only the linear portion 
is utilized. MAZE [ REf. 8 |] computer program is used to generate specimen 
model, grid, define boundary conditions, and also create the input file that 
can be immediately use for NIKE2D program. ORION [Ref. 9] is the post 
processing of NIKEZD which can be used to create contour curves of stress, 
displacement, etc. The observation of uniform stress in each element will 
be made. The observation of using half of the entire model by symmetrical 
will be made, which bases on the idea that the model can be divided into 
smaller elements. when compare to the entire model. The linear 
relationship between external loads and the local stress in each element 
will be verified. The plate without notch is used as bench mark to compare 
the reliability results to the plate with notch. The specimen model 
dimension bases on specimen dimension which using in the shear test 
experiment. The shear model that involves with the multiaxial stress 
which can not yet be treated in this thesis. Future work requires 
generalization of current results to multiaxial stress. When such extension 


is made, realistic prediction of composite structures will then be possible. 


H. STATISTICAL MODEL OPERATION 


A. REVIEW OF WEIBULL WEAKEST LINK FAILURE THEORY 


consider an infinitesimal volume element 6V at point x and subjected 
to stress state Ojj. The probability of failure &f of volume element 6V 


is assumed to be given by: 


6f = ¥ oj (x) ] 8V ( 2.1) 


wher e: Y = failure parameter which js a function of 9j which is a 


function of position, x. 

It is assumed there are only two possible states for the volume 
element, &V , either the volume element.survives the loads or fails. From 
this assumption the probability for survival or reliability is simply: 

BR=1- bf = 1- HL oj lx] 18V (2a) 

The basic assumption In weakest link theory is: If any volume, 6V, in 
the total structure volume, V , fails, then the entire structure occupying 
volume V fails. 


Consider the structural volume, VV, to be divided into N volume 


elements §V, , with each of the volume elements located at point xy. 


Then the total volume is given by: 


where: Ni = 00 
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The probability for the volume V_ to survive, R, under weakest link 


theory, is given by: 


R=lim T SR, =lim 1 {1- ¥f oj le,] 1 8Vy } (2.4 ) 
where: k=1,2...N, and 
N+ oo 


Equation 2.4 can be equivalently written as: 


InR= lim Z in{l - VE Gj [xy] $V, } (25 ) 
where: kK=1,2...N, and 
N= 00 


The quantity ¥.8V, is small and from the Taylor series expansion 


InQi+x)~x for x <<] 


Equation 2.5 becomes: 


InR=lim = -{ ¥f Gj jlxy] ]8Vy } | (2.6 ) 
or R= exp -( fy ¥ ola) ]4V } | (27) 
where: K=1,2...N, and 

Ni -» 00 


Equation 2.7 is the basic equation derived under the assumption of 


weakest link theory. Note that the stress distribution need not to be 
uniform; the material can be generally anisotropic and the material can be 
non-homogeneous. The equation is limited in that it describes only the 
reliability of the entire structure and provides only a quantitative 


description of where the structure will fail. The reliability becomes: 
R = exp - { );/8 ie 
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where: « = shape parameter, and 
8 = scale parameter 
In the next section, basic theory will be used to obtain the 
formulation of size-effect for one dimensional, non-homogeneous, 


uniaxial stress. 
B. SIZE-EFFECT FOR NON-HOMOGENEOUS STRESS, ONE DIMENSIONAL, 
UNIAXIAL STRESS 


The value Of oy, and Bi are estimated by the experimental 
procedure from the given m numbers of sample length L!. For a sample 


length L , the evaluation of o and 8, can be implemented according to 


the following procedure. 
There are j_ links in length L (Figure 2.1!-a) each links are not 


necessarily of equal length, but the stress in each link is homogeneous, ie 


the stress oh within the p"" link satisfies the following condition: 


xO? - (fi) KO,AKI/I, S$ € CZ 
where: € Can be any arbitrarily small value 
ly = length of p!” tink 
p=) ee 
The reliability for each link is: 
Rp(o,?) = 1 - fo?) = exp -(,P/B,)% (29) 
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Figure 2.1 (a) One Dimensional Links Configuration 
(b} Two Dimensional Specimen Configuration 
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There are (In/1) unit links within Ip , the reliability of each unit link 


within Ip is: 
Rp(o,P) = exp -(o,P/8 )™ (2.10 ) 
where: « and 8 for unit link are material constants 


The reliability for each link Ip is: 
Rp(oy?) = TT {exp -(o,P/8) ry 7M = 1, 2.2... (aD 
= exp (-(I,/IMo LS 
where: (I,/1) is the number of unit link in |, 

The reliability for then entire length L Is: 

Ry = Ry (oy ) Ro(G5) Rz(03) Ryg(og) ..- R;(o)) 
= exp {-(Ih/1(0 1/B)°F exp {-(12/1)(o 2/8). . «exp {C1 /ICo 4/8) 
= exp ~{(1/1)(o 1/7 B)™ + (12/1) 2/B)™ +... + (i /1)(o 5/8) 


Re = exp -£ E [(Ip/1(o F/B) } = I, acs) “Ca 


Under uniform stress condition within an element, let o, Is the 


p 


element stress which is constant within ie length element and Oy iS 


calculable from a given boundary condition (B.C.) characterizable by {P}, 
(eq, a load vector). The stress due to the same B.C. with a different 


magnitude P is: 
OpIP] = op {PH (PAP) ( 2.12 ) 


where: P can be of any arbitrary magnitude in linear system, or 


P is of neighboring magnitude of local linear system. 
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op[{P Hl is expressible from structural mechanics calculations, 
(eg, a finite element method.) 
In order to reduce the value of o,P which is a function of P in 


equation 2.11 to a single random variable {P} , the substitution of o, PIP] 


by oO, PI{P}] is made within the same manner as the expression in equation 
2.12. Thus, the equation 2.11 becomes: 

RL = exp -( E[(Ip/o yPUPYIPYXP/B)I)} = sph... 

Expressing the above equation into the Standard Weibull form: 

Riiexp -{(P/ BeAr } 

Then: oq = 

B, = B/ (Z[p/IMo PUPI/PH/™ 5 p=t2... | 

Let: = F=(Z1(I,/1) (6, PLPHI/(P)IE plas, 5 = aia 

Thus: By = BF 

That is, given a structure of length L subjected to boundary condition 
parameter {P} ( which may be of a single or multiple load Pi, P2,... 
vectors characterizable by a single parameter {P} ) which give rise to a 


unidirectional stress 0, computable from the methods of mechanics ( eg, 


finite element analysis ). I are element lengths which are seqmented as 
small aS necessary to assume the stress to be uniform. Upon sampling n 


structures measuring P,, failure loads. 
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oy , £, are shape and scale parameters ( in dimension of load 


estimated from DIES and structural size L. 


o« , § are shape and scale parameters ( in dimension of stress which 
are geometric independent material parameters ). 

F is a conversion factor calculable for structural ( finite element ) 
analysis. 


For special case, homogeneous, One dimensional uniaxial stress: 
From: F=(E((Ip/1) “*(o,PUPHAPII IY sp=12...) 
For o,?=o and {P}=P_: 
a Oe a 
={(0/P,) {2 (Ip/1) 4% 37! a) ae 
= (0/P,) | (E/1) 
Thus: 
R, = exp -{ [ (P/B)/(a/P,) | AL/1) “1% J} 
= exp -{ (L/1) [ (o/B)(P/P,,) 1™ } 


This recovers the well known classical result demonstrates the 
consistency of the current formulation is true. 

For each given non-homogeneous boundary condition, a corresponding 
correction factor F can be calculated by finite element such that the size 
effect can be predicted. By weakest link in length X:; it is inferred that 


the plane of facture is perpendicular to xX. That means the failure plane 
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has no thickness. Since the one dimension specimen has no width, fracture 


width is not a parameter. 


C. SIZE-EFFECT FOR NON-HOMOGENEOQUS STRESS, TWO DIMENSIONAL, 
UNIAXIAL STRESS 


This is the extension of one dimensional, uniaxial stress to two 
dimensional uniaxial stress. The definition that uniaxial stress implies 
failure has zero dimension in X direction is adopted for one dimensional, 
uniaxial stress , which has been described in the previous section. For two 
dimensional, uniaxial stress, failure automatically requires definition of 
the failure plane because the specimen has dimension in Y direction 
(Figure 2.1-b), and it is no longer conceptually a point as in the one 
dimensional case. In the similarly procedure as in one dimensional case, 
the specimen is divided into sufficiently small areas within which the 
stress in direction of force is assumed to be uniform. It is seen that not 
only the length of the specimen is divided into small element lengths but 
also the width must be divided into small element widths in order to 


provide element areas which are satisfies uniform stress assumption. Thus: 


KO yP - 1p f xO, dx/1, < €, for all |, 
where: eal 2. ny 

tae wal YFxdY/Wa S & for all Wa 
where: Gulez ...K 


Weakest link in X direction is defined by Eq. 2.8 and unit length in xX 
direction is defined by Eq. 2.9. 


Ay 


Weakest link in Y direction is: 
Rqlo,P] = exp -£0,.°9/8, 4} pq 
where: o,P9 is stress in the gih element, and q=1,2...k 
This means that if any one of the git fail, all elements fail. However, 


the unit (metric) width can not be arbitrarily defined. It is visualized that 


within each metric area, there is one domineering crack (largest crack 
perpendicular to o,). When subjected to o,P4, the domineering crack 
extends without bound. From fracture mechanics: 

ko =O a, 


where: k . is a material constant and can be independently measured. 


Cc 


do is weakest link in Y direction, which can be used as the unit 
metric width. 
Thus: a,P49 = (k-/0,P9)4 for element pq! 
As well as the division er length element Ip Into unit length, the 


division of width element is necessary. In stead of unit length 1, a.?9 is 


used in order to indicate the critical width in each element. thus, the 


reliability of entire element is: 


Rl? 4 = exp -{ (Ing/ I(w nq/ ach 1o,PA/ Ree 
= xP -£ (Ing/IIW pa/{ke/FyP%)21(0,,PV/ 8) } 


= 2x0 -£ (Ipq/IN(W pq/kc2(o,PN2°%/ BX } 
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where: pq =1, 2... jk ( subscript is changed from two dimensional 
array into one dimensional array.) 
The reliability of entire specimen (which is the structure) is: 
Re = Ry Ro Rs .. . Rix 
= exp -{ E (Ing/I(W 5 g/ KANO, Pq)Ztaxy gx } pg=1,2... jk 
From the linear relationship between external load and local stresses: 


0, PAP] = 0, PAP) (P/{P}) 


wher e: oPSI{P}} is element stress which is calculated under certain 


loading boundary condition while the magnitude of load is {P}. 


o,P4[P] is element stress which subjected to the same 


boundary condition as above but different in magnitude, the magnitude is P. 


Thus: Rp =exp -{£ (Ing/ WOW g/ /k_2)((P/{P} })0,,PSE{P}] ZX Bx | 


Pq * 
where: pq =1,2... jk 


Express it in term of the Standard Weibull for the structure, thus: 


Re = exp -{ (P/ Be) GzZ5)) 
where: oe x + 2 


Be = (PIB6O/%*2) (5 I G/M pq/ke2Mo,PAEP IN? JOM Coer2)) 


where: Bae 12. 2 uk (Za) 

The values of « and 8 of material are altered by the boundary condition 
and the geometry of the specimen (or structure) changed but they do not 
depend on the changed of the magnitude of the load. Therefore, the shape 


parameter, « and the scale parameter, B can be solved for any given 


OS, 


structural geometry and boundary condition. The finite element method can 

be used to calculate the stress distribution of the specimen by dividing the 

specimen into small elements within which the stress element is uniform. 
For the simple limiting 2-D case, the specimen without notch, the 


stress distribution for entire specimen is uniform. Thus: 


0,P4 = 0, = constant 
Be = ((P}/o4) Blo/ot 2h, (2/42) { 511 Wo gl Mt) 


= {P} k(2/%G) glor/oE) ig, all/oe) 


where: A = EU IpgWpq | ,pq=i,2... jk 
XE = «+2 
Thus: Ae = exp -{ (A/k,) (P/{P}) (,/ 8% }°E ( 2.15 ) 
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Il. APPROACH TO THE PROBLEM 


A. SPECIMEN MODEL CONFIGURATION AND COORDINATE SYSTEM 

The specimen model has 1.99 ” in length and 0.968” in width is align 
in the y-z plane, (cartesian right hand coordinate system is used). The 
configuration and the coordinate of the specimen without notched and the 
plate with 90 degree notch configuration are shown in Figure 3.1. The 


fiber direction which makes the orthotropic material having large elastic 


modulus ( E,, ) is considered to be along the y-axis (Figure 3.2). 


B. FINITE ELEMENT PROCEDURE 
|. Material Constants 
The finite element program requires defining the material 
constants [ Ref. 7 ] for orthotropic elastic to be used in stiffness matrix 
calculation. The material constants of typical graphite epoxy composite 


which used in this thesis are: 


E,, = 18.65 Msi (Mega Pound per square inch) 
vz) = 0.0339 
039 = 0.4 


2) 
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Figure 3.1 Dimension of Specimen Models without Notch and with 
Notch 
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Specimen Models Coordinate System 
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Gio = 0.84 Msi 


2. Boundary Conditions 
The specimen model is subjected to the displacement boundary 


condition (B.C.) in order to be able to use the entire specimen model in the 


calculation. An original problem the specimen model is subjected to the 
displacement B.C. value wu, (Figure 3.3) on the right hand side in the 
positive y-direction and the left hand side is subjected to the displacement 
B.C. value u, In opposite direction. In order to fix the specimen model 


from rigid body displacement, the left side of the model must be constraint 


in the direction of load. Therefore, the displacement value on the right 


side of the model becomes two time of The The transform of the 


displacement B.C. is shown in Figure 3.4. 
3. Data Planning and Purposes 
NIKE2D is two dimensional axis symmetry finite element code 
that the thickness is assumed to be infinite. The plane strain analysis is 
used by assuming no displacement along x-axis and no rotation around y and 
Z axis, which it is satisfy for this purpose. There are eight finite element 


program data output to be discussed in this thesis: 


3.1 Specimen model without notch subjected to the 
displacement B.C. value 0.01” 


3.2 Same specimen as 3.] but the displacement B.C. value is 
0.001". The purpose of these two cases is to demonstrate that there is no 
stress concentration along displacement direction and the linear 
relationship between displacement and local stress. 
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Figure 3.3 Displacement Boundary Condition 
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3.3. Specimen model with notch subjected to the displacement 
B.C. value 0.01” and the stress calculation bases on four-point integration 
(Ref. 7) in order to observe the uniform stress in each element. 


3.4 The same as 3.3 but the calculation bases on one-point 
Integration at the centroid of each element. 


3.5 Half of the specimen model with notch subjected to the 
displacement B.C. value 0.01” and four-point integration is used. This case 
will be used to compare with the entire specimen model in order to insure 
that the half specimen can be used because of the symmetrical 
configuration. 


3.6 Same as 3.5, but one-point integration is used. 
3.7 Same as 3.6, but the displacement B.C. value is 0.004. 


3.8 Same as 3.6, but the displacement B.C. value is 0.001. 
One-point integration results will be used to be input data for 


Statistical program. 


C. STATISTICAL PROCEDURE 

After obtaining the results from the finite element program, 
four-point integration results have to be considered to insure that the 
tensile stress in each element is uniform or almost uniform. If they are 
not uniform, the size of the element need to be reduced. However, due to 
computer memory limitation, absolute uniform stress may not be 
attainable. One should realize that what spatial domain that the element 
size should be reduced and what area the element need not be changed. The 
symmetrical of the model can be considered to be used only half of the 
entire model in order to increase the number of elements. 


oy 


The statistical size effect in two dimensional, non-homogeneous, 
uniaxial stress program in fortran 77 is written to be used to calculate the 
reliability of the specimen models by using the finite element results and 
shape parameter, « and scale parameter, 8 from several numbers of certain 
Sample size. There is One program including: 1) the specimen model 
without notch, 2) the entire specimen model with notch, and 3) the half 
(top side) specimen model with notch. | 

This program is presented in Appendix A. There are 12 output data 


files which are shown in Table 4.3. 
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IV. DISCUSSION AND CONCLUSION 

A. EXAMINATION FOR UNIFORMITY OF UNIFORM STRESS WITHIN 

ELEMENT 

The four-point integration are obtained in entire and half specimen 
models subjected to the displacement boundary condition (B.C.) value 0.01. 
For the full specimen model, the element number 544 (shade in Figure 4.1) 
is Considered, the tensile stress values can be read in figure 4.2 ,where the 
Stresses of this element are tabulated under rows 2173, 2174, 2175, and 
2176 (relatable to the element by dividing the last number by 4 recovering 
544) and the tensile stresses are 0.03663 Msi, 0.03738 Msi, 0.01933 Msi, 
and 0.01851 Msi (Figure 4.2) respectively, which demonstrated that the 
stress within element is not exactly uniform but the differences are small. 
The element number of the half specimen model at the same position as 
element number 544 of the full specimen is 442 (shade in Figure 4.3) and 
the tensile stress can be read in Figure 4.4, where the row numbers are 
1765, 1766, 1767, and 1768 and the tensile stresses are 0.03667 Msi, 
0.03742 Msi, 0.01937 Msi, and 0.0186! Msi respectively, again they are not 
exactly uniform but the differences are small. The element number 192 in 
Figure 4.5 (full specimen model) is outside the notch region, the tensile 
Stress values can be read in Figure 4.6, where the row numbers are 765, 
766, 767, and 768 and the tensile stresses are 0.0689] Msi, 0.06888 Msi, 
0.06852 Msi, and 0.06855 Msi respectively , which shown that the stress 
within element is not uniform but the differences are more uniform than 


element in the notch region. The element number of the half specimen 
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model at the same position as element number 192 of the full specimen is 
230 (Figure 4.7) and the tensile stress can be read in Figure 4.8, where 
the row numbers are 917, 918, 919, and 920 and the tensile stresses are 
0.06906 Msi, 0.06905 Msi, 0.06871 Msi, and 0.06872 Msi, while they are not 
exactly uniform but the differences are also smaller than the element in the 
notch region. In principle the stress in each element may approach 
uniformity by decreasing the size of the element. This will increase the 
number of the elements and ts limited only by the computer memory and 
computation time. The observation of symmetrical model is obvious from 


the stress values of full specimen model and half specimen model. 


B. OBSERVATION OF LINEAR RELATIONSHIP BETWEEN EXTERNAL LOAD 
AND LOCAL TENSILE STRESS » 


The tensile stress of the element number 460 (Figure 4.3) of half 
specimen model which subjected to the displacement B.C. value 0.00) ts 
0.002475 Msi (Figure 4.9) and the tensile stress of the same element 
which subjected to the displacement B.C. value 0.004 is 0.0112 Msi (Figure 

4.10). The ratio of tensile stress is almost the same as the ratio of 
displacement, the different value bases on the non-uniform of the stress 
around the notch region and can be accepted. When the element number 230 
(Figure 4.7) of the same model, which is far from the notched region the 
ratio of the stresses (Figure 4.11, displacement B.C. value 0.001 and 
Figure 4.12, displacement B.C. value 0.004) and the ratio of the 
displacement are almost the same and the difference is small compare to 
the different value for the element in the notch region. Therefore, the 


linear relationship between external load and local stress can be accepted 
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whether the element is in notch region or not. Since the element is small 
enough to be satisfy the uniform stress within element, then the uniformity 


and the linearity can be attained. 


C. STATISTICAL SIZE EFFECT RESULTS DISCUSSION 

Equation 2.13 and 2.14, calculated the reliability of the specimen 
model based on statistical size effect in two dimensional, non- 
homogeneous, uniaxial stress. For application in these equation, the 
displacement boundary condition need to be converted to the average stress 
along the boundary. The definition of the average stress in two dimension 
is: 

(P= (2 LW qo,PI }/ (2 wWoq) 
where: pq (element number) = 186, 187, 188, 189, 190, 200, 205, 210, — 
215, 220, 225, and 230 for half specimen model (Figure 4.7), and 

pq = 163, 164, 165, 166, 167, 168, 171, 174, 177, 180, 183, 186, 
189, 192, 465, 468, 471, 474, 477, 480, 483, and 486 for entire specimen 
(Fiqure 4.5 andFigure 4.13). 

The values of stress in each element are shown In Table 4.1 and Table 
4.2 respectively. 


The statistical strength of a structure is characterized by the scale 


parameter Be in equation 2.13 ( same as BETAA in computer program, 
Appendix A). A change of the §¢ value is therefore related to the strength 
change of the structure. The values of Br (equation 2.14) for plate with 


notch, where the values of «, 8, and kK. are 10, O11 Msi, and 0.00! 
4] 


0.001 respectively are shown in Table 4. 3-A. Table 4. 3-A also shows 
that the variation of the values of Bre which are calculated from different | 
values of external load (PBAR) from different values of displacement 
boundary condition (u’). The values of B- which calculated from half 
specimen model seem to be higher than Be value from the full specimen 


model. The reason is that the number of elements for full specimen is less 
than the number of elements for half specimen model, therefore, the 
accuracy for full specimen model could be lower than that for half 


specimen model. For the purpose of reliability calculation for the 
structure, one can average the Br values from one set of finite element 
data output from different values of load under the same boundary 


condition, therefore,for instance the average value of Be is approximately 


0.030. Table 4. 3-B shows the values of 8, where the values of «, and k. 
are changed to 20, and 0.0001 respectively. The Be Values are reduced but 
the difference of the average 8, values is not large, because of the 
compensation of higher « and lower k, values. This means that the 


changes of « andk, values do not affect much to the $- value but « affects 


the reliability of the structure (equation 2.13). All the reliability 


calculation results for different values of Br ,where the specimen models 


are Operated in the same loading range are shown in Appendix B. 


A2 


In the case of plate without notch that shown in Table 4. 3-C and 


Table 4. 3-D, the values of Br from different values of displacement 


boundary condition are exactly the same. The reason is that the uniformity 


and the linearity affect the accuracy of the calculations. The value of 8 


from the plate without notch (0.0439) is higher than from the plate with 
notch (0.028) under the same boundary condition (u’ = 0.01), which is a 
reasonable result. 


In addition, the results show that the material scale parameter, 8 is 
converted to Be by the geometry of the structure and the loading boundary 


condition, but it does not depend on magnitude of the load. The results of 
the specimen model without notch confirmed that the uniform stress 


assumption Is valid and known result can be recovered. Between results of 


notched specimen model we showed that 8 remains constant for different 


load magnitude because the uniformity and the linearity conditions are 
attained. As expected we confirmed the notched model is less reliability 
than the plain model under the same loading condition. Given available 
experimental data of a composite material, we can evaluate the scale 
parameter 8B and shape parameter «. Using these results we demonstrated 
that the procedure of analyzing the effect of a non-uniform state of stress 
obtained from the finite element method can be post processed by 
Statistical reliability formulation accounted for size effect, resulting on 


the prediction of the reliability of the structure. 
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DSF = 2. t3GE-=G8 
TIME: @. Q@GE +a 


TEP Te 


Ca he Poa hat 

PAAR ee _fe aa f 

SRE SES of op sbee, Assh sh a po 
aja 243 24 5 













2G5 9cec S24 247 268 26 9 240 ro age eq3 cas mS é e 2 247 P27 B 2 249 24a 
249 nate é92 e493 294 25 Seas 247 28 ete 268 e 24 1 24 2 24 ho 


243 294 2395 236 24? 238 2 p24 CAD 2A. She 249 2445255 246 247 246 
8 ayo 2; 2 edi ja e 2 ed3 B24 4 2a5 ¢d6 2 ’ 248 cagee 2: 2 pa: 242 


Figure 4.1 Element Numbers of Entire Specimen Model at Notch 
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ELEMENT NO. L 


w 


SIGY 


S1G2 


a1G2 


SI1G6Y2 


2159 9.49124375 9.0124375 YJILYLIDEH91 MV, 5795E MNS ALI 245E M02 HN, 997JEHD2 
2150 9,9124375 0.0124375 Jied1ANE MDL -0,9U58E “03 OL1096E 902 917,92 30E=02 
215] 0.9124375 0.9129375 I.5028E-01 O0.5551E-0% O,2212E"02 20.9784E=02 
Ploe 0.0124375 0.0124375 9.5051E=-01 9,.9SIUE#03 0.,.2e3588FE "02 =-9.1059E=9]1 
2153 0.9124375 0.0129375 YI.371HESN1 HO0.1K4O3E S03 0,1295E 202 -0.1056E-01 
2154 0.0124375 0.9124375 YJI.3H91E $91 99,5359F 913 O0L1118F 202 “0.9745 E=02 
P1165 0.9124375 4.0124375 0.524SE<-03 9,1N9UES02 N.2307E “902 -0.1032E=-01 
Clbo 0,7124375 0.0124375 0.5279F-01 O0-1638E=G2 O.2554EH0e eDOLIIITE=01 
2157 0.0124375 0.0124375 O.3218E901 OLY77REHDNZ OLI3S9LTEH02 -0,1109FE 0} 
2168 0.01245375 0.0124375 I W3183E901 “0, 7994EMRN4 NLI14UE M02 ©0.1022E=-01 
2159 39.9124375 0.0124375 O.4oOudEW01 OL1757E#02) 0.2395E"02 -0.1072E=01 
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4.12 Stress Distribution Values of Half Specimen Model Based 


on {-Point Integration, where u’=0.004", at Right 
Boundary Region 
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Figure 4.13 Element Numbers of Entire Specimen Model at 
Bottom-Right Boundary Region 


56 


TABLE 4.1 


AVERAGE STRESS FOR HALF SPECIMEN MODEL WITH NOTCH 


NUMBER 


186 
187 
188 
189 
190 
199 
200 
205 
210 
219 
~—220 
oa 


{ 


ELEMENT 


0.057000 
0.057000 
0.057000 
0.057000 
0.057000 
0.024875 
0.024875 
0.024875 
0.024875 
0.024875 
0.024875 
0.024875 
0.024879 





0.009071 
0.009358 
0.009496 
0.009994 
0.009750 
0.008758 
0.00849) 
0.008179 
0.007840 
0.007500 
0.007195 
0.006968 
0.006895 


0.036280 
0.037420 
DOE TIS) 0, 
0.038190 

0.038280 
0.035030 
0.033960 
0.032710 

0.031350 

0.029980 
0.028760 
0.027840 
0.027880 


WIDTH, IN. | U’ = 0.001” | U’=0.004” | U’ = 0.01" 


0.090670 
0.090350 
0.094830 
0.099370 
0.099970 
0.087550 
0.084860 
0.081700 

0.078260 
0.074810 

0.071720 

0.069400 
0.068220 


= PBAR" 0.008717 0.034850 0.087045 


Average boundary stress 


oy 


TABLE 4. 2 


AVERAGE STRESS FOR FULL SPECIMEN MODEL WITH NOTCH 


| ELEMENT 


Oo 


NUMBER 





WIDTH, IN. 





163, 166 0.095000 
164, 167 0.095000 
165, 168 0.095000 
171, 465 0.024875 
174, 468 0.024875 
177, 471 0.024875 
180, 474 0.024875 
183, 477 0.024875 
186, 480 0.024875 
189, 483 0.024875 
192, 486 0.024875 


y Msi. 


U’ = 0.01 


0.091560 
0.094940 
0.095670 
0.087600 
0.084860 
0.081660 
0.078170 
0.074670 
0.071540 
0.069200 
0.068030 


PBAR 0.078025 
*% 


Average boundary stress 
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TABLE 4.3-A 
STATISTICAL DATA SUMMARY 


FOR PLATE WITH NOTCH, «=10, B= 0.1, k.= 0.001 


MODEL HALF . FULL 
oe | u'=0.001" | U'=0.004" | U'=0.01" | U’=0.01" 
PBAR. 0.008717. | 0.034850 | 0.087045 | 0.087025 | 


0.030982 0.030749 0.031110 0.028880 














TABLE 4. 3-B 
STATISTICAL DATA SUMMARY 


FOR PLATE WITH NOTCH, «= 20, B=0:1, k.= 0.0001 


u’=0.001" | U’=0.004" U’=0.01" 
0.008717 0.087045 | 0.087025 


Average boundary stress 



























0.029122 0.029721 | 0.028826 
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TABLE 42 5-C 
STATISTICAL DATA SUMMARY 


FOR PLATE WITHOUT NOTCH, «=10, B= 0.1, k.= 0.001 


PBAR 0.009512 0.094910 


o> Fo 0.043948 0.043948 










TABLE 4. 3-D 
STATISTICAL DATA SUMMARY 


FOR PLATE WITHOUT NOTCH, «= 20, 8 =0.1, k,= 0.0001 


PBAR- 0.009512 0.0949] 


a er ror 


Average boundary stress 
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APPENDIX A 


STATISTICAL STRUCTURAL RELIABILITY PROGRAM 


Creme heme RORP Rare meme Eee REM ORR ehhh ROR eR 


Ce 


Ceere rater are keeraekrekindkeee metre h ah tithhaetkekaeeeenakkreneeke ee 


E 


Cre Re RRR RARER REET RERERERERRARRARRERRARARARAA ERR EREREN 


Oo 


C 
G 
C 
c 
C 
S 
C 
C 
C 
C 
G 
G 
C 
C 
C 
G 
C 
ec 
C 
C 
C 
C 
C 
C 
C 
C 
e 
C 
G 
c 
c 
C 
C 
c 
e 
c 
s 
C 
c 
C 
C 
C 
C 
C 
c 
C 
c 
C 
c 
C 


PROGRAM: RELIABILITY CALCULATION OF PLATE MODEL BASED ON 


LINEARIZE SIZE EFFECT, NONUNIFORM UNIAXIAL STRESS. 


RELCABILITY FORMULATION IS EXPRESSED IN THE 


" STANDARD WE[BULL FORM “, 


DATA INPUT: QUTPUT OF FINITE ELEMENT COOE NAME * NIKE20" 


VARLABLE: KC = 


BETA 


ARakkeRKR ERK RRR aR RRRE RARER RE RR eH ARH ARERR RRAAKRAHARARKRARAAA 


MATERIAL CONSTANT FOR PLATE THAT IS 
SUBJECTED TO UNIFORM TENSILE STRESS. IN 
THIS CASE, ASSUME IN EACH ELEMENT THE 
LOCAL TENSILE STRESS 1S UNIFORM, ~ 

= SCALE PARAMETER FROM EXPERIMENTAL RESULT 


FROM SEVERAL NUMRERS QF SAMPLE. 


ALPHA = SHAPE PARAMETER FROM EXPERTMENTAL RESULT 


PBAR 


FROM SEVERAL NUMBERS OF SAMPLE. 

= AVERAGE STRESS ALONG BOUDARY SIDE OF THE 

Plate THal SUBJECTED TO THE EXTERNAL FORCES. 

IT MUST BE CALCULATED SEPARATELY AFTER OBTAINING 
THE STRESS RESULTS FROM "NIKE2D" PROGRAM, 

NUMBER OF ELEMENTS OBTAINING FROM MESH GENERATDR 
PROGRAM NAME “MAZE” WHICH IS THE PROGRAM THAT 
CREATE AN INPUT FILE 70 BE USED IN “NIKE2D" 
PROGRAM. ALSO BOUNDARY ELEMENT NUMBERS CAN BE 
CREATED FROM "MAZE" PROGRAM. 


DOUBLE PRECISION F,FA,FF,.P(1000),RAC1000),BETAA 
DIMENSTON XL (2368), XWC2368).S1GY (2368) 


REAL KC 
BYTE AC30) 
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aAaMAmMANNnNnAaNNM 


1 


mw 
wm 7 


23 
24 


TYPE*,/ENTER FILE CONTAINING STRESSES DATA INPUT’ 
READ(5,5) A 

FORMAT(30A1L) 

OPEN(UNIT=4, FILE=A, STATUS='OLD') 

OPEN(UNIT=6, FILE='SS.DAT', STATUS='NEW') 

WRITE (6,14) 

FORMAT(3X,'DATA OUTPUT FOR LINEARIZE SIZE EFFECT 


tUNIAXTAL TENSILE STRESS') 


TYPE*, ENTER {1 FOR PLATE WITH NOTCH 


he OTHERWISE FOR PLATE WITHOUT NOTCH! 


READ*«,NC 
IF(NC.EQ.4) THEN 
WRITE (6,2) 

ELSE 
WRITE (6,22) 
END IF 
FORMAT(20X¥,°FOR PLATE WITH NOTCH’) 
FORMAT(48¥,"FOR PLATE WITHOUT NOTCH’) 
TYPE*, ENTER THE VALUE OF DISPLACENENT BOUNDARY CONDITION’ 
READ*«,¥Y 
WRITE(6,3) U 
TYPE*,*ENTER 4 FOR ENTIRE SPECIMEN INPUT 


te OTHERWISE FOR HALF SPECIMEN INPUT 


REAOD*-NS 
[F(NS.EQ.1) THEN 
WRITE(6,23) 
ELSE 
WRITE (6,24) 
END IF : 
FORMAT ({5xX,'ENTIRE SPECIMEN MODEL TNPUT'?) 
FORMAT(16X,'HALF SPECIMEN MODEL INPUT') 
FORMAT( 7X, 'WHERE THE DISPLACEMENT BOUNDARY CONDITION [S ',F8.4) 
TYPE*, ENTER THE KC VALUE’ 
READ*+KC 
TYPE*x,' "ENTER THE BETA VALUE' 
READ*,-BETA 
TYPE*«,'ENTER THE ALPHA VALUE' 
READ*,ALPHA 
ALPHA2=ALPHA+2.0 
PFzO 35: R 
TYPE*, ‘ENTER THE PBAR VALUE’ 
READ*,PBAR 
TYPEx, 'ENTER THE NUMBER OF ELEMENTS' 
READ« -NE 


REAO ELEMENT LENGTH, WIDTH, ANO TENSILE STRESS FROM INPUT See 
CALCULATE BETAA WHICH [S BETA THAT CHANGED WITH GEOMETRY OF THE 
SPECIMEN MQDEL AND CERTAIN LOADING BOUNDARY CONOITION BUI IT LS 
INDEPENDENT QF THE VALU= OF LO4D UNDER THE SAME BOUNDARY CONDITION, 
DO 200 1=1.NE 

READC4,20) XLCI)/XWCI).SIGYCTI) 


IF(SIGY(I).GT.0.0) THEN 
F=(XL(L)*xWOL) (KC €*#2.0) 


se *(SLIGY CTI) s*AL PHA2 ) 


FFF Rt 
ELSE 
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Amr 


109 
20 
2s 
30 


FF=FF 

END JF 
CONTINUE 
LF(NS.EQ.1) THEN 

FF =FF 
ELSE 

FF=FFa2.0 
END IF 
FA=FF ee (-1./ALPHRO) *PBAR 
BETAA=BETA**(ALPHA/ALPHA2) *FA 
WRITE(6,10) KC 
FORMAT(8x,'KC =',F 13.8) 
WRITEC6,11) BETA 
FORMAT(6X,'BETA =',F13.8) 
WRITE(6,12) ALPHA 
FORMATC5X,'ALPHA =',F13.8) 
WRITE(6,13) PBAR 
FOR“AF(GxX,*PBAR =',F13.9) 
WRITE(6,14) NE 
FORMAF(8X.,'NE =',[4) 
WRITE(6,15) ALPHA? 
FORMAT (Qx,'ALPHA2 =',F13.8) 
WRITE(6,25) BETAA 
TyP&x, ENTER THE FIRST EXTERNAL LQAD' 
READ#,P(4) 
PINC=P(1)/30. 


RELIABILITY CALCULATION BASED ON LINEARTZE SIZE EFFECT. 


DQ 100 j=1,100 
RACJ) =DEXPC-(P (J) /BETAA) «*AL PHA?) 
WRITE(6,30) J,P(J),-RACI) 
P(J+1)2P(J)+PINC 
CONTINUE 
FORMAT (8X,.2F11.0,E12.90)  . 
FORMAT( 5X, 'BETAA =',F43.8) 
FORMAT(9xX.‘°S =',34,3X%.'P =',F13.8-3X.'RA =',F 335.8) 
STO? 
ENO 
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APPENDIX B 


DATA OUTPUT FOR LINEARIZE SIZE EFFECT UNIAXIAL TENSILE STRESS 


DATA OUTPUT FOR LINEARIZE SIZE EFFECT YNIAXIAL TENSILE STRESS 
FOR PLATE WITH NOTCH 
WHERE THE DISPLACEMENT BOUNDARY CONDITION LS 0.9654 
HALF SPECIMEN MODEL INPUT 


KC = 0.00040000 
BETA = 0.10090000 
ALPHA = 20.00000000 
PBAR = 0.900871700 
NE = 482 
ALPHA2 = 22.00000000 
BETAA = 0.92912212 
J = l P = 9.90871799 RA = 1.90009090V009 
Js 2 Pos 0.90990757 RA = 1.99990900 
Jos 3 a 0.90929913 Mone 1.09900009 
J= | P= 0.90956979 Ra 2 1299000000 
Ji 5 Ps 0.9N987927 Ras 1.00000000 
J = 6 Ps 0.9101639483 RA = 1.00000009 
J = 7 Ps 0.91046049 RA = 1.90000N000 
Jos A Pp = 9.91979997 RA = 1.000099000 
Js 9 Ps a) Vl RA = 1.990909090 
J = 1U Ps 0.113321" RA = 1.90090009 
J 21 Ps 9.91152267 RA s 1.09000000 
J == Le P = 0.91191 323 RA s 1.09900009 
J = 13 P = 9.01220 589 RA =s 1.009000000 
224 Ps 0.912939%37 QA = 0.99999999 
J = 15 es VeAIETAYDS RA 0.39999999 
JI2 16 a= 0150755" an) = 9.99999993 
tos? icy Ps 0.014$34507 RA = 0.999993996 
J = 18 P = 9.013555603 RA = 0.99999994 
gos A519 Ps J.91394720 QA = 9.9999999]1 ~ 
Js 20 P 5 9.01423777 24 2 0,99999985 
Jz ?1 Pos 0.91452933 RA = 0.993999977 
=e P= P.,o14atagny RA = 0.999299965 
SPS 23 y-s JLOVSLO947 oes 9.9997999ud6 
Js e4 P os Vet 54903 aA = 0.99999918 
joes 25 Pr os 9.915959959 ao. = )},9F999RT77 
J = 26 po OL9ISIALI7 RA = 0.99999R815 
Js 27 Ps Vole fies Ss aA 0.99999725 © 
eS 5 lee G a= 9.91630230 Bas 0.99999595 
J = 29 Ps 0.915935 287 PA = 0.9999941)6 
j 7s. 30 2s VI Pia sas = 0.939999135 
8 P os vie saa 2A i 9.9999A747 
paz ase 2 = 991772157 RA = 0.99998198 
Js 323 P 3s O59 td0 Ts Be is 0.99997423 
Js 3a a9 = V.013830570 QA = 9.99994330 
ee ae Ps 09.01359527 RA = 0.9999uead19 
J = $6 p= Ne US2R S45 QA = 0.99992712 
Js ee ps jj, 097410 QA = 0.99989A04 
J = 38 P= J.919%6 797 ies 0.99985SAU045 
J = 39 Pos 0.01975353 QA s 0.99920337 
J 2 «0 Ps 2292994919 24 = 0.99972A90 
Jz atl y= 9.02033767 QA s 0.999452797 
J = ae os 9.92993 1)23 24 s 0.99939175 
JI = 4&3 Ps 9.0299792930 RA s 0.°99930879 
J = 44 Ps U2 ee les 7 or a 0.999637 3 
Js 4§ Ps 0.02150193 94 3 0.99873722 
J = 46 ef 9.79179255) wa = 0.9933037A 
J s9a7 Pos 9.0229H397 RO = 0.93977306) 
J 2 48g Ps %002237303 RA = 0 .99697554 
Je es = 1202266420 RA = 0.99592476 
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DATA OUTPUT FOR LINEARLZE SIZE EFFECT UNIAXIAL TENSILE STRESS 
FOR PLATE WITH NOTCH 
WHERE THE DISPLACEMENT BOUNDARY CONDITION 18 0.01400 
; HALF SPECIMEN MODEL INPUT 


KGn> 0.00010000 
BETA = 0.190090000 
ALPHA = 20.00009000 
PBAR = 0.08704500 
NE = 482 
ALPHA2 = 22.000900000 
BETAA = Q0.02972093 
Jos l Ps 0.)09871790 RA = 1.00000909 
J = Pag Ps 0.99900757 RA =s 1.99000000 
J = 3 Ps 029992984 3 RA = 1.09000900 
Jos a) Ps 9.0995A870 RA = 1.009909009 
ci Sa S Ps }.09937927 RA = 1.0290090000 
Jos re) Ps 0.91910983 RA = 1.000000900 
Jo= if P = 0.919046040 RA = 1-09000000 
Jos 3 a ae 0.J1975:97 PL s 1.90V000000 
ee 9 = PSO TI OATS 3 2A = 1.00000000 
J = LU Ps Veott $3210 RA = 1.90000009 
Js it as O.91152207 RA = 1.09000090 
J = ie Ps 0.911913823 RA = 1.00000000 
J= P = 0.01220359 RA = 1.00000000 
J = Ts P= 0.012949437 RA = 0.99999999 
Jos os 0.01274349%3 RA = 0.99999999 
J = 16 Ps 0.91397350 wo = 0.99999999 
= ae aa 0.01336507 RA = 0.99999998 
J = P = 9.01355563 RA = 0.99999996 
Ios 19 a 0.01 594720 AA = 9.99999994 
j= 20 So es C02 5777 RA = 0.992999991 
Js 21 o = Ceti 5 295 5 RA = 0.99999985 
w= ee Ps 9.91491994 RA = 0.99999974 
J = 25 Ps eI91S1994N7 QA = 0.99999996 
= 24 a oe 0.015 79003 RA = 0.99999Q94A 
J 2 25 Ps 2.9159996) RA = 9.99999921 
i= fe6 2 = Q2015 78117 ras 9.99999HR2 
Jos) ce P= eee 7 17 3 FA s 0.929999825 
Js 29 Ps 99165023 QA s 0.99999741 
JI = 29 a — V.U16°5247 RA = 0.99999629 
ee 30 os O.O1714 343 Pa 49.99999447 
JI= $31 Ps Y,O174954NN RPA = 0.99999200 
gee, Se Ps Se 7 ie A 7 RA = 0.99998a8au8 
J = .33 P= F.01891513 RA 0.99998353 
Js 34 PS 0.01830570 RA = 0.99997658 
G7= 35 pees, 9291859527 RA = 0.99996689 
Nees. 6 P= DVI AA5B3 QA = 0.99995 343 
io. oF a= ere 7 te PA = 9.9999 3484 
J = 38 Ps V.019%0797 RA = 0.99990929 
J = 39 Ps 0201975353 RA s 0.999A7T4U34 
j= 40 P= 0.029)4910 RA = 9.99982675 
J= atl Ps 9.072033947 Anes 0.99976224 
Jj = 4e Ps Hse s se 4 aA = 0.999675)]A 
Js 63 2) 20902092940 RA s 0.99955817 
ios 4a P= 9 svelcil lst RA = 0.99940158 
Js 65s Pps 9.)2150193 YA = 1.99919284 
Jo= U6 Ps 0.02179250 RA = 0.9949157N 
ay 7 Ps Dates as 07 2A = 0.995549t4b 
J = 48 Ps VeMAAST Fy 3 RAs 0.998N 6624) 
Jos 4g Pos V.0225n4An wA s 0.997%3222 
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DATA OUTPUT FOR LINEARIZE SIZE EFFECT UNIAXIAL TENSILE STRESS 
FOR PLATE WITH NOTCh 
WRERE THE DISPLACEMENT BOUNDARY CONDITION [5 0.090040 
HALF SPECIMEN MODEL INPUT 


KC = 0.000190090 
BETA = 0.10000000 
ALPHA = 20.00000000 
PBAR = 0.034850900 
NE = 462 
ALPHA2 = 22.00000000 
BETAA = 0.02929771 
J = 1 = Vous 1700 RA = 1.00090000 
Jos Pd Ps 0.96900757 RA = L.90000000 
J = 3 Pos 0.099729313 RA = 1.09000000 
i 4 Ps 0.9995 397) RA = 1-00000000 
wos 5 P = 0.09947927 Se ae 1.99000000 
Js 6 JS 0.u1016785 RA s 1200000000 
J = 7 y = 06971946940 RPA = 1.00000000 
Jos A Ps 92910975097 RA = 1.900000900 
eas 9 Pos 72911040155 QA = 1-090900900 
Js 19 cae 2911353219 Ra = 1.09000000 
Jos. Vial Ps J201152267 RA = 1.09000000 
jz sive se ae Peel oises RA = 1.-00900000 
Jos YS a= DeJ12293An RA = 1.00000000 
2) a 14 es 9.01249237 as 0.99999999 
J = eas P= JeM12S749495 RA = 0.93999999 
Jo2> 16 » = 0.91397550 Rh = 0.99999998 
Js lee Pe s 0.01336507 RA = 0.99999997 
J = ws pos ye Soo o0 5 24 = 6.999239995 
JI=s 19 ps od1 $97 SNH PA s 0.993999992 
JI= 290 Ps IN 257 7 on S 9.99999987 
0 eee P= Ve Olea e 8 55 ahs = 0.9999998(0) 
j=, £2 as de1491789N RA = 1,.992999969 
Jey Gs pos PN15199387 ea Ss 0.93999953 
Jee ees Pos 1,015 Fons Ra = 0.999990]98 
eS 2S a 0.915592 61 2a 6S 0.99999K92 
J = 26 — QeAMtSIeil 7 tA = 0.999998 34 
Gus V27 Pos )<) tee 7 lize RA = ).299999759 
Sys 32s Ps 9.714552 350 Ras 0.99999645 
J = (238 Ps VeJl6452 37 RA = 0 ,93999479 
J = 3 P =z e017 125s aA os V,99999240 
ip aS 31 - s Ue017 234 PA = 0.9292998993 
J-= $e Ps O28 Li i23a7 RA = 0.9999H42] 
Js 335 oS = OQ.O1891513 RPA = 0,9999774? 
Jos tq Ps 0.914303 70 RA: s 0.939967909 
Js 35 Ps 0.01859527 (Qa s 01.93999S"6uU 
J = gs Paes Pe VI5%H583 FA s 9.99993615 
Js 37 P 3s Q.91917740 PA = 0.99991067 
J = 48 Ps 0.9194797 QA s 0.93947564 
j= 52 P= 0.91975355 eo. 2 0.99982772 
Jo 40 p= 1.020048719 RA = 0.9997624A 
‘ae = 9292983967 PA s 0.999657405 
J = Ps 0.02093 )7 $3 RA s 0.99955470 
= een Ps 0.329920 40 RA = 0.99939431 
6 ee ee Ps ).92121137 RA = 0.99917957 
Ji 45s Ps 2202156195 QA s 0.99989357 
Jes te os 9.92179259 VA = 9.99851375 
J = o7 P = V2 e sn RPA = O,99801137 
J = &8 » = OS ee 5? aos xA = 0.99734976 
Jj = 49 ps FeDFAA5H 4A Da fs I.9354A1 50 
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DATA OUTPUT FOR LINEARIZE SIZE EFFECT UNIAXTAL TENSILF STRESS 
FOR PLATE WITHOHT NOTCH 
WHERE THE DISPLACEMENT BOUNDARY. CONDITION [5S 0.0010 
ENTIRE SPECIMEN MODEL INPUT 


KC = 0.00010000 
BETA = 0.10960000 
ALPHA = 20.00000009 
PBAR = 0.00951200 
NE = 200 
ALPHA2 = 22.90000000 
BETAA = 0.05180000 
Jos 1 Ps 9.00AR717900 RA s 1.9090000000 
j= 4 ge s 0.90909757 RA = 1.099000009 
Js 3 Ps Jen 99298 5 QA = 1.00909990 
j= | = 0.9995R49A70 RA = 1.99000000 
J ee S Pas 9.999487927 9A s 1.900900000 
ae 6 = VINLV1IO IR 3 RA = 1.09009N000 
J = T Ps 2.910469 40 RA =s 1.090090900 
J os R Ps ..9L975997 RA =s 1.00000000 
Jos 9 Ps AS04153 QA = 1.090000600 
Jos 10 Ps 0.U1133210 RA =s 1.09090900 
j= Ps 09.91152267 2A sc 1.000000900 
JI= 12 Ps ve 3s1s25 RA s 1.90000000 
Js 4S os 0.91220389 Sp 1.00000V00 
J os 14 Ps 9.012998 37 RA = 12909000000 
JS LS os I.IDI2TANDS RA = 1.053900 9 
Js Vie Ps 9.01307550 RA s 1.900090000 
- Ve eos 157135436607 RA &S 12.00000000 
Je 13 Ps 0291355563 RA s 1.09900000 
Jos 1 a= 0.01394729 RA ss 1.09009000900 
Jos 29 a et A257 7 7 RA s 1.909990000 
J =) 2 2S 9.014952933 oe 1.09099000 
J = 22 a ea O00 RA = 1.00000000 
= 28 o = Devt S1 9907 ho = 1.9990000Q0 
fs Wea = 00154923003 RG = 1.09900000 
j= 225 a= V.9b5 590 Ht RA = 1. 000%0000 
J = 26 os (eo rat i rA = 190000009 
JI = e@7 oh Weoloe? 17 § PA s 130700000 
JS = 28 Ps uv .U1630239 == 1.06009000 
Js 29 Ps 01945297 PA s 1eH90nyH0N 
J-= 939 p= DedL714 343 RA = 1 .000090909 
Js eo Ps VOL 79 3400 KA = 1.00Q000U009 
J = ‘32 Pos Mele te ts wA s 1.09000NQ9N 
Js £53 Ps 0.01491513 Ra 1.9000 99000 
Js 34 Ps 0.912303579 RA = 1.900009000 
J) .2 233 P= ea a9 S27 Qa = 1.99009000 
Js 50 Pos Yo. 1339535 PA = 1.09N009000 
J = 30 Ps O.9F917740 QA s 1.900N0000 
J = 38 a= 0.71945797 PA = 1.9009)00 
J = 39 ys Nee? 3355 RA s 1.09000000 
J = 40 Pos AN2098919 ey = 1.0909000900 
y= el aes 0.729353767 RA s 1.09000000 
J = te Bs Mege 1557025 ai c 1.90900099 
J = U3 o = PeO20972Z IA es 1.09090N00 
J = Qa Ps Hee 21 2 bs 7 RO s 1.909000900 
J = a5 aa Gevetsyt ys vA = 1.9900N0090 
J = U6 P= 9.9217929% RA = 0.99999999 
J =u Pos 0.2P2N KENT QA = 0.99299999 
i= SS pars 1.02237 4563 PA = (6992999949 
Js 49 pos V972PSR ION RA Ss 0.999999990 
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DATA OUTPUT FOR LINEARIZE SIZE EFFECT UNIAXJAL TENSILE STRESS 
FOR PLATE WITHOUT NOTCA 
WHERE THE DISPLACEMENT BOUNDARY CONDITION IS 6.0190 
ENTIRE SPECIMEN MODEL INPUT 


KC = 0.00010000 

BETA = 0.10000000 

ALPHA = 20.000000006 

PBAR = 0.09491000 

NE = 200 

ALPHAZ = 22.009000500 

SETAA = 0.05180009 
J = 1 Pie VONKBTITYO RA = 1.999000000 
J =z 2 P = 9.00956757 QA =z 1.990900N00 
‘i 3 P = 0.:0929913 AA = 1.NNNDDHNND 
J =z 4 P= OL9N95A979 QA = 1-090000900 
iS 5 P = 0.9098 739e7 RA = 1.99090000 
J = 6 P = O.91TO 1694S RA = 1.990090000 
J = 7 P= G6.010446950 aas 1.09000000 
es a Ps 0.01)075097 RA = 1.0909009690 
J = 9 P = vert rotts3 Pa =z 1.69)09000N 
Jo2= 10 P= VoOL133210 RA = 1.000000900 
JS Stal P = FSG 1S2 267 QA = 1.99000990 
J= te Ps 0.011791 323 PA = 1.909000009 
J = 4G P = 9.91220 540 RA = 1.990000909 
J = Bie as Oe 1124e1 57 QA xs 1.99060000 
Js Bie Ps 001278993 202 1.99000090 
J = Gis Pos 9.01307359 QL = 1299000000 
J] = “ae pos 0913505597 PA = 1.9900000UN 
J = 18 Ps 9.91355593 PA = 1.900090009 
Jj os 19 — 0.,01394720 24 2 1.070n00N0 
J = Bo P= p eal ee ee ar 9 24 2 1.UIJ09GN00N 
= 21 Ps 2.914952 7535 aa = 1eGONN00000 
J = Bee P= eJ1 IIL AID P4 = 1.0C000009 
)2 23 P= OSULS1TO747 WA = 1.990)90000 
Js 24 P = 001540995 9A = 1.990000U0 
ios ee P= Je) 959559 ees 1.090900000 
jes "S26 a= (Poots aie 7 Po 3s 1.99090090000 
jes 327 Ps Vounlo2I{73 vA = 1.0900 YUNGD 
Js 28 $$-P = 0.163623) RA = 1.99900000 
as 29 ae YodtHag247 RA = 1.9000690000 
A ee) Pp = e917 eS a 1.29000%u0 
J = $1 Pp = Letts 9 14 yA t.9900G0080 
J= $2 Bee Peli letoe On = 1.0200NG00 
j= Gees re N9.0PAN1TS13 2A 2 1.9390900090 
J = 34 Pos Sirs soa RFA = 1.09990909 
J 3 SS a cs 0.91359527 he = 1.09900N00 
Js 36 P= VoVPARR S43 RA 1.90009000 
= 3s a= OeVLIIL77 40 PA = 1.09000000 
J = 38 p= VoIpata7a7 Qi = teganoronnn 
j= foe = 7.01975453 Ra = 1209900900 
al uf es 2e VENNYQIL!) PA s 1.VINNNVDN 
y= ay P= MP0 33I957 ie 1.999000000 
Js ae os 9209295328 RA = teedoguanunan 
Jos &3 P= 2eI2092 AY 74 = Leva O0000U 
4 ee ud > = 2 92 121-1 cane ce 1.99090003) 
Jo= as Ps ‘Veith 2 1 50s: ay = 1eN90900U0)) 
i = US ae= Vous | aes wA = 0,9a39999999 
V2 5am Pp = 3}. 2 27A3-)7 RL = 0.97929999 
J oz 4g oe =: Nee 257535 Or = 1.99999999 
|} = 2a9 -  s Veovaeon ten is = ..99999994 
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DATA OUTPUT FOR LINEARIZE SIZE EFFECT UNIAXIAL TENSILE STRESS 
FOR PLATE WITH NOTCH 
WHERE THE DISPLACEMENT BOUNDARY CONDITION IS 
ENTTRE SPECIMEN MODEL INPUT 


0.0100 


KC = ¢.00010000 
BETA = 0.1900090000 
ALPHA = 20.00000000 
PBAR = 0.08702500 
-NE = $92 

ALPHA2 = 22.900000060 
BETAA = 0.00288262 

J 3 t eS 0.908717590 FA = 1.990N00000 

J = 2 eec= 9.09900757 RA = 1.009000000 

J = 3 P = 0.909929913 QA s 1-09009N0C00 

jz " Pos 1.99958979 RA = 1.00000000 

j= 5 P = 0.9N9473927 RA = 1.000909000 

J = 6 Ps 0.01901463983 i= 1.90090000 

Joz 7 ps 0.01046940 RA = 1.00090N00 

J = Q Pos 92919750997 RA = 1.09000900 

Jos 9 Ps 9.01104153 RA = 1.N9099000N 

Jos 10 Ps 9.911332190 RA = 1.90000000 

Js 11 Pos V.91162247 RA = 1.00000000 

J = ile ae = 9.01191 323 RA = 1,.0000)900 

Jos 1$ Pp = 1.01229390 RA = 0.99999999. 

Jo = ie Ps A eo Sus 7 PA = 0.99929999 

J 2 €5 Ps Ve91 278193 RA = ..99999998 

2 = "16 Ps Mot 507 550 Qi = 0.939999997 

d= P = A213 30507 RA = 0.99999995 

J = 8 Pos 0.91355543 ens 0.99926992 

J= 19 — 0.0139472N RA = 0.999999RAR 

J = 2% 2 = orld es7 77 Rh = 0.999999R1 

J = 2! Ps mur 1523.55 AA = 0.9999997] 

15s. S22 PoP = OSDTIAEIOG 94 s 0.9999995% 

J -3) hes P= Ces) Fa 7 PA = 2.99999930 

J = 24 Pos eos 190 5 vA = J.99999994 

Jos 325 coc =: NS71559 760) OA s 0.999998u40 

Jos 26 o = Neos 4 117 WA s 0.9999976% 

J = 227 Pas A.V1H27175 oy ee 0.99999643 

j = oA ps 0.N10356239 RA s 9.92999474 

Js 29 Pp = I.ulodS24%7 PA = 1.99999228 

Js Sv Pos Cae? 1G 303 RA = 0.9999BA76 

J = ds Q.017134090 RA = 0 .99998373 

J = 3e Ais Veo 7257 RA = 1.99997bb60 

io eas P = 0.91KO1513 RA = 0.99996653 

J= 34 2) RS 0.918309579 QA s 0.9232995241 

t= 83> P= Qala 52 | RA = 0.9999327{ 

J = 30 Pos 0.90183 9AR543 IA sc 1.999990536 

J 2 7 P = OOIDtL77 40 RA = ..99986759 

Js ese Ps 0.01940/797 RA = 0.999R1567 

J s 3a Pas Det 5753 rate 0.99974544H5 

J7= 43 os 0.0203491u Po = 0.99954795 

fee es 1.9°2034495?7 QA = 0.9995 lonAs 

j = We P= 0202053923 PA = 0.999540 01 

jz 4% ps 8.92N92SI9G eis 0.99910233 

j = ue P= VoN2141137 RA = 0,798 7AN2P7 

J = ies P= PeIAIS0193 HA = 0.999360 39 

J = 46 Ps O02 hi 2250 YR = 0,99779774 

J = -e7 Qs Hee eCasy FA = VL99TNS3R3 

J 2-246 p= Weare Sees wi = P.FFOUTUIOR 

) = ag OF 2c Devecwns | a ie W QQ VTRATO 

P= (S50 = os Wee es 2 lee 4 CAFS, any 
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DATA OUTPUT FOR LTNEARIZE SIZE EFPECT UNIAXLAL TENSILE STRESS 
FOR PLATE WITH NOTCH 
WHERE THE DISPLACEMENT BOUNDARY CONOZTION [5 0.0010 
HALF SPECTMEN MODEL INPUT 


KG.e= 9.900100000 
BETA = 6.1909900000 
ALPRKA = 10.00000000 
PBAR = 0.00871700 
NE = 482 
ALPHA] = {(2.00000000 
BETAA = 0.03058189 
Jos l Ps O.9CATIT00 RA = 0.9999997] 
J = 2 o> 3 9.009f0757 RA = 0.99999957 
je 3 no PeQNI29R{3 aA = 0.9999993A 
J = 4 Pos 1 .0095AR970 RAS 0.999999190 
Sites 5 Ps GO, 0nss7 927 RPA = 0.99999A871 
Js 6 2 = -0.01016983 RA =< 0.99999817 
J = 7 P = 9.91%%6040 RA = 0.99999744 
Jp = a pecs On 59 7 Sf = 0.99999644 
Josz 9 Ps VOLINGIS3 RA = 0.293999509 
J= 19 eo = Senet. 55.2 10 RA = 0.99999330 
Jos 11 P= 0.01152267 RA = 0.99999092 
J = le P= 9291191323 RA = 0.99998779 
Js 13 o s 9.91220389 RAs 0.99998369 
} =a Ps 9.01219437 RA = 0.99997837 
J eral p=: 9.9127A393 RA = 0.99997150 
JS = Jo pie 0.01307550 PO = 0.99996268 
J =e P = 9.91330507 PA = 0.999951 42 
J = ele Ps 0.013554503 RA = 0.99993712 
ee =) ys ),9'394720 RA = 0.92991904 
| 20 , s YoU) 425777 Pa ss 0.99989632 
Jos 1 oa VLA1LN5 2233 PA s 0.9°9986787 
J = 2e P= O.91391899 RA = 9.999R3044 
Js 233 Fos FQISLIIAT 9A 5 0.9997RRU7 
J = 234 oa 20915990903 cs 0.99973415 
J => P= 0.01S59060 QA = 9.99966732 
J S924 Ps Ds oP ely Wa = 9.99958539 
J fee Ps Oa 2a le § RA = 9.999%HS34 
J = 28 a= 0.916356239 Pao = 9.99936359 
J= 29 Pos 0.916%5247 aa = 0.99921595 
Js 30 p= OC91L714343 th = v.9990375S2 
jGa5 St Ps Veen 2A = 0.99832257 
Jee ise = 001 Pes 7 aA = 0.998560443 
Jj = B33 Ps VOLH01513 RA = 0.99825539 
Jos 34 P = \ UVES o 579 PA = 0.997FABH49 
J iS Se Pos Y.01859527 RA = 09.99744740 
Js 36 do s Y.J1AZA543 RA = 0.99592624 
J 25s P = OSOLIL7ZT NWO 2 aa 0.995309 56 
J = 38 Ps 21.919 lo797 Q4 = 0.99558119 
OS Sys pae= 9.91979453 RA s 0.99472362 
Jos a0 o s 9.929049190 RA = 0.99371655 
Js” 4 4.920383 967 a= 0.99253675 
J = ue a= To 2NS Sie 5 RA = 0.99115800 
ee Ps Pe DMAVIAIRY 24 2 0.9895567 
= es i 0.021 24+137 wA = 05939764135 -- 
J = 45 P = 9.92150193 9a = 0.9959551251 
Js ub P s D027 Via) RAs 0.98300219 
Jos a? Y= O.0IDIWASN7 7A = 0.99019314 
J = ae pos 0.0223874303 RA Ss 0.97576349 
abe 4g - = YCMIPSH IAN iS V.97TAIPSIS 
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DATA QUTPUT FOR LINEARIZE SIZE EFFECT UNTAXIAL TENSILE STRESS 
FOR PLATE WITH NOTCH 
WHERE THE DISPLACEMENT BOUNDARY CONDITION IS 
, HALP SPECIMEN MODEL LNPUT 


0.0040 


KC = 0.001900000 

BETA = 6.19000000 

ALPHA = 106.00000000 

PBAR = 0.03485000 

NE = 492 

ALPHA2 = 12.090000000 

BETAA = 0.03074947 
Jos 1 Ps 0.90871790 RA = 0.99999973 
J os 2 P= OS 009007 S37 QA s 9.99999960 
a 5 ct Me NIV e2 OR ts RA s 0.999999%2 
Js 4 Pos 9.999583 7(0) RA = 099999915 
J == 5 P= 0.90997927 RA s 2.29999H79 
jis 6 P= YQeUVLOLS9B3 QA = 0.99999829 
a 7 ea X.01040040 RA = 0.99999760 
los g 8 s 91.91075097 RA s 0.99999666 
J= 9 P = Oot lo4tss aoe 0,99999549 
Js 10 Ps C2011 35 200 RA = 0.99999372 
Jos 11 P 3s 9.01152257 9A = 0.99999150 
Js 12 Ps Yo911791323 RA = 9.99998856 
JS ies Ps 0.01220 380 GA = 0.99998473 
j= Ps 0.01249437 RPA = 0.99997975 
j= oS Ps 0.012738193 RA = 0.99997331 
Jos 16 Ps OSU Ss 07 530) 24 = 0.939996505 
J = “Wa ps }.01336507 ws 9.99995450 
J = 18 P= Veo 555 5 ans 0.9999411] 
jz WW Ps 0.9139472u RA 0.99992%14 
J = 20 Ps .9.,01423777 WA os 0.999990290 
Js 21 Pos 0.914523233 QA = 0.99997626 
Js 22 P= V.0TaEs1 890 RA s 9.9FI9DR43ZN7 
J = 23 —— O.01S1U947 RA = 9.99980190 
Jo= 2% Ps UVeN1S49NI03 Gas 0.99975193 
J s 25 a= 9.01559909 2A = 0.9996ARY SZ 
J = an Ps Veo Ui 2A = 9097961179 
des. ae Pp Peete 7175 RA s 0.99951 R00 
Jos 28 ps 329103023) QA 02 .99940397 
J = 29 P = DeUte3%5237 RA = 0.99926579 
= se ys AUL7TLE3ZI3 244 2 2.2990985A 
J = 31 a= F-91795 45990 ca ).99B89726 
J = Se Ps ea eae Ares) 74 RA = 0.99955549 
Js 35 Ps Ue 8 a a a= 0.998356U 
ts ee ee 0201539579 Rai = 0.9930 2N051 
y= 5 oe = 229018592627 R4 = 0.99760)923 
J => 36 Ps Ree eed 2% adi = 9.997121N6 
J = (oy Ps Y.91997739 ao = 0.99554321 
J = 38 Ps QVs91%416/797 ro. = 0.995361 01 
d= 739 Oe 02791975353 2A = 0.99595770 
J = 0 ee O.02004910 QA = 0.739G011420 
J 2 at ps 1602033967 Pa = V.993ZNN0R8ARO 
J ee Qo = 092053923 QA O.99L71 690 
ts pa pes 9.22092 989 HA s VIINA1NG7 
J = a4 P = Dethe 1 eit 47 Ya s ).98845873 
J = es 28S VIPAT SINI93 RA = 0.995425B1 
J = % —— Jee 7 254) QA s OS IetGre 55 
J 2 67 o = Vob2ecAR su OA = er es Se 5 
Js 466 aa Je92257305 IA = V.9T9?P2P 21 
J) = 49 P= W,92A5H 12% QA s O.971462159 
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[ 
DATA OUTPUT FOR LINEARIZE SIZE EFFECT UNIAXIAL TENSILE STRESS 
FOR PLATE WITH NOTCH 
WHERE THE DISPLACEMENT BOUNQARY CONDITION IS 0.0100 
HALF SPECIMEN MODEL INPUT 


KC = 0.990100000 
BETA = 0.10000000 
ALPHA = 140.00000000 
PBAR = 0.08704500 
NE = 482 
ALPHA2 = 12.00000000 
BETAA = 0.03111023 
Jos 1 P = Q.008371790 RA =s 9.99999977 
J = Po Ps 9.90900757 Qa os 0.99999955 
Js 3 PS 49,99929413 RA = 0.99999909 
J = a Ps 0.089584870 RA =s 0.99999927 
Jos 5 P = 9,00997927 RA = 0.99999A95 
Jos b Ps VeILVIMIAS RA = 0.99999851 
J = 7 Ps N.919%6940 RA = 0.99999791 
J = a P= 0.901075997 Ras 0.99999710 
J os 9 Ps Q9.VIt9G1s$ RA = 0.99999691 
J = 10 BS 0.9113321°) RA = 0.99999454 
Jos Hal Pos 0.011592267 aA> = 0.99999261 
J- =) We P= OeOLLI1323 RA = 0.99999006 
J os 1 Ps 0.91220389 RA = 0.99998672 
a] So ee P= 9.01249437 no = 0.99998239 - 
J™= = to P = O.0L27AR9S vA = 9.°99997680 
J = lo P= 0.V1307550 RA = 0699996962 
J 22? RPo= 0.013%36497 RA = 0.99996045 
J = 18 p= 91.01365563 RA = 0.,99994AH0 
J>= 19 Ps 0.0139472% aA = 0.9299 34008 
J => 620 P = 0.91423 707 PA = 0.999915S5A 
Jo= 21 p= 9.91452933 RA = 0.99999 242 - 
J 2 2? &Pese 0.901431999 RA = 0.99986357 
Js 825 Pos 0.915190947 RA s 0.99932777 
jz 24 P= 91.0154h403 RA 0.9997R354 
its 25 Ps 0.9155900 Pa = 0.99972912 
J = 26 oD = 0.01598 17 AA = 0.999564 1 
J = 27 Ps Y.01677175 A, = 0.9995A8094 
ios 25 Ps 0.01556230 RA = 0.99948180 
Is, o20 aS 0.916952487 RO = 0.99936158 
ioe=e aoe! Ps VOT 7 (aes 9A = 0.99921627 
Js $i eas 0.91723400 RA = 0.999041 22 
Js) 3e P= 0.91772457 RA s 0.99823109 
ea P = 0.01901513 RA = 0.99857930 
J = ja Ps 0.01830570 RA = 0.99827883 
I=. rod Ps 0.01859527 —- RA = 0.99792117 
J = 36 P=: O.N1LAAIAKK3 RA = 0.997%9601 
Jos 37 Ps 0.91917740 RA = 0.99699403 
J = 38 ea = 0.91946797 R4 = 0.99640063 
j= 3a Ps v.01975353 RA = 0.99570183 
Jo= a9 Ps 9,020N48I71N RA = 0.994R809R 
A ame) a= 9.92933907 RA = 0.99391915 
J = &e Ps 0.920533923 RA = 0.992T79%86 
Jos 43 2 = 22092092980 RA = 0.9919838') 
J = &4a BPs - 0.021 21837 RAs Os98S9905R55 - -- - = 
j= wus P = 9.02150193 WA = 0-IFALBBLCY 
P= 46 Ps 0.92179250 RA = 0.9961382) 
i= a7 2 BS 0.92204307 aa s 0.99376971 
) = as Ps 0.922%7303 AA s 0.98103944 
Js ag a 02922564279 Qaos 0.9778994NS5 
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DATA OUTPUT FOR LINEARIZE SIZE EFFECT UNIAXIAL TENSILE STRESS 
FOR PLATE WITH NOTCH 
WHERE THE DISPLACEMENT BOUNDARY CONOTTION IS 0.0100 
ENTIRE SPECIMEN MODEL INPUT 


020197 Sits 3 RA 
9.020904919 RA 


O.99955721 
0.9875 4645 


KC = 0.00100000 
BETA = Q.19000000 
ALPHA = 10.00000000 
PBAR = 0.98702590 
NE = 592 
ALPHA2 = 12.90000009 
BETAA = 0.02887978 
oe t P s 0.90871790 RA = 0.99999943 
Jos 2 p= 0.09990757 RA = 0.9999991I5S 
i oes 3 Ps 0.09929451 3 oy 9.9999948745 
4 es u Ps 0.)0958879 RA = 0.99999821 
Jj = S P= 0.900987927 RA = 0.99999743 
Jos 6 Ps 0.01016983 RA = 0.99999636 
Jos 7 P = 0.01946I40 oe 9,.99999490 
J = 4 Ps O,01) 75097 RA = 0.99999292 
J = Q Ps OL OULI4S1S3 RA = 0.99999025 
Jos 10 Ps 9.01133219 RPA = 0.99998456028 
J = Bia P = O20 11 s22o7 RA = 0.99998195 
J = Jj2 Ps 0.0t191323 RA = 0.99997572 
J os 13 er = 0.01229380 RA = 0.99996758 
Le 14 P= 0.01239%37 RA = 0.99995799 
i ee 15 eae 0.91278493 R4 s 0.99994334 
J = ts Pcs 0.913907350 RA = 0.99992581 
Jo= 17 ears 0.01336507 2A Ss 0.99990342 
Js VA P= 0.01355503 RA = 0.99987498 
J 2 19 Ps 0.91394729 RA = 0.99983905 
Js 20 Ps 0.01423777 RA = 0.99979387 
Jos tf Ps 00145233 5 nes 0.99973733 
Js 22 Bs 0.01491390 PA 2 0.99960089 
Js 23 P= 0.91510947 ie TE 9.9995799 
v= a7 = Hgts 10003 RO = 0.99947153 
J = 28 oa = 0.9155906" RA = 0.9993 3849 
j= 26 P= 0.91598117 oes 9.99917S87 
2 ae 27 Baas 0.01627173 Pa = 0.998977N4 
Js 2 ea 0.9160350239 RA = 0.99873513 
J = 29 P os 0.01635287 Ra = 0.998441 A} 
jis 3a P= 201714454 3 Sane 0.998038737 
Jf = $1 Ps 0.01733400 "A Ss 9,.997660047 
ues 132 P = 0.017724%57 RA = 0.99714793 
Js | 35 P = O.01LaN1 513 QA = 0.99653447 
Jos 33 Ps 0.913830579 RA = 0.99580?745 
Jo= 35 Ps 9.91959527 ‘RA = 0.99493150 
JI = 36 OC s NLOLAXR SRG PA = 9.99389AR25 
es ees o = N50 19 1 7949 RA = 0.99257590 
JI 2 38 Ps 0.019%6797 RA = 0.991233453 
= Oe P= = 
J = 2 Ee = 
J = @i P= 9.02933907 R4 5 0.°%93S2167P 
J = Ste P= 0.U2)5302 8 aS 0.98 249774 
Tove: as el 0.92092980 Sa a 0.97933e61 
J os gu- Ps 0.02121-137 Ra = 0.97565798 
ure aS P = V.02150193 pa = 2.971403515S 
J = t6 y = 0.02179259 RA = 0.75548977 
Jo 2a Pics 9292204307 RA = 0.95993139 
J = Un Ps 7209237363 Pa 9.95433317 
J = 49a os Vet22 sou OS 9.94589176 
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DATA OUTPUT FOR LINEARIZE SIZE EFFECT UNIAXTAL TENSILE STRESS 
FOR PLATE WITHOUT NOTCH 
WHERE THE DISPLACEMENT BOUNDARY CONDITION IS @.0010 
ENTIRE SPECIMEN MODEL INPUT 


KC = 0.00100000 
BETA = 0.10000000 
ALPHA = 19.00000000 
PBAR = Q.90951200 
NE = 200 
ALPHA2 = 12.00000000 
BETAA = 0.943948035 
Jp | Ps 9.998717 390 RA = 1.09000UU00 
jJ os 2 P= 0.90900757 RA = 0.99999999 
J = 5 a= 92999294813 BA = 0.939999999 
J = 4 os 06.99095A970 RA = 0.92999999 
J = 5 Ps 0.009987927 Ra = 0.9999999A 
J = 6 Ps 0.91010783 RA = 0.99999998 
Jz 7 o> = O.N1LU45040 Rh Ss 0.99999997 
Jos @ aes 0.01075997 9A = 0.99999995 
Jos 9 Ps Cont 0a $ QA s 0.99999994 
Js (td Ps Os01153210 no = 0.9999999] 
Jo= It Ps 0.91192267 RA = 0.99999988 
J = Se Ps 0.91191323 pac Ss 0.999992984 
jos 13 Ps N.VL22N 389 RA 6S 0.99999979 
J = Wel oO = 0.129937 Ri = 2.99999972 
J = Ge Pos 9.012738493 RA = 2299999963 
JI => {06 vos 0.0130755Su RA = 2.99999952 
J = Ay a= 0.013346597 RA = 0.99999937 
J = 90s os 0.01355543 QA = 0.93999919 
Jos 19 Pos 9.01394726 Ria s U.99999896 
J = geo ES 0.014257 77 WA s 9.999998 665 
| = een Fs Jeo 5 as 5 RA = 9.99999830 
J =aeze 9 3 9.014391 990 QA = 0.99999784 
Js 23 os 9.91510947 QA = 0.99999727 
Js 24 o = Y91546993 P4 s 9.99999657 
J z= a5 28 $e 0155 706") aa = 0.9 es 57 | 
Jz 2» 88 VeUISPALIZ = FA =) 0,99999465 
Js 27 es J. Vlb2 tae 3 2A = 0.999993 36 
Js 325 Ps 0.01m30239 pA = 9.999991 79 
Jo 29 pos 0.01045297 QA = 9.99998989 
J = 30 oc = yl 7 1s 3 Qa = 0.99998759 
J = Pp = JeIL 7839900 RO = 9.99998481 
J = 382 Ps 9501/7 72557 vA = 0.999981 4A 
J = 5S P= 0.91891515 Sia = 0.99997749 
J >= m4 P= 0.018 3Ge7" an i= 0.99997273 
Jo= 35 Ps V.918S59527 RA = 0.99996705 
J = es oP = 0.01 2832R583 Rios 0.99996031 
Jose 37 > = 0.919177490 RA = 9.999952 $4 
J= 38 a 1.019%0797 QA 0.939994291 
Jo= 39 Ps 9.01975453 RA = 0.99993180 
J = 40 P= 9.0299499710 QA = 0.99991 875 
J = 4@t P= 1292933967 2A = 0.99990 343 
J = 4y2 els 1. Ue vs Se s aoe .9F9RBSS 1 
Js 43 Ps 0692092989 24s 9.99986159 
a = Ga Pos V OF, 41+37 RA 2 Y.IIIRUDN2Z2 
Jos 45 pos Jev213I193 RA s 0.99981189 
J = 4h 2 ES VIL 7 9254 aa. = 0.9997791)1 
J= 47 o = 02298347 rao 0.999743095 
} = WS Pos 1.2731 5 24 2 999959696 
Jo= a9 P= 1692250425 eA = 9.°9954622 
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DATA OUTPUT FOR LINEARIZE SIZE EPFECT UNIAXIAL TENSILE STRESS 
FoR PLATE WITHOUT NOTCH 
WHERE THE OISPLACEMENT BOUNDARY CONDITION IS Q.0100 
ENTIRE SPECIMEN MODEL INPUT 


KC = 9.00400000 
BETA = 9.109000000 
ALPHA = 10.09000000 
PBAR = 0.09491000 
NE = 200 
ALPHA2 = {2.90000000 
BETAA = 0.04394902 
Js 1 p= O.908717900 ea, = 1.9N000000 
J = 2 Ps Ge 700757 24 3 0.92999999 
Jos BS Ps VeUUIZIK] 5 Qa 2 0.99999999 
Jigs i Ps 0.9095827)) RA = 0.99999999 
J = 5 o = 0.90997927 as 0.99999998 
Jos 6 Pos Je9L91H 933 RA = 0.99999998 
Joz 7 Ps 9291090456940 es 0.99999997 
Jj a Paes ‘AO SULO7SA97 PA = 0,999999095 
a 9 P os O.OLL 24153 oe = 0.99999994 
J o<= 1Q Ps Oe tT 353520 RA = 0.9999999] 
a tt fs 0.01152257 RA = 9.9999998A 
J = Sie P = 0.01191323 QA s 0.99999984 
J = at Ps MeOte2 3550 RA = 9.99999979 
Js 14 P = Pete 194650 | 2A 0.99999972 
j= 1S eis 9.01279 493 vA s 0.932999%963 
J = 16 Ps Oe 30.7550 RA = 0.999999Se2 
J = ee eos 9.913565%7 RA s 0.99999957 
J = 1A Pos Veul $55505 RA = 0.92999949 
jos 19 Ps QO.0139572°9 RA = 0.999799 89% 
J = 29 Ps: geola23s7 77 po 0.939998654 
j=" ot Pos 0.314523 53 RA = 9.99999839 
J = 22 Ps Ore) 31200 PA = 0.999997R4 
J = 23 — PeOLSt O97 PA s 939997 27 
hos, §24 Bos 9.01539903 Qh =: 3.99999657 
Jos 25 os Jeo 3s 7 15°) RA = 9.99999571 
= 26 — 0.019 394117 a. = 0.997999 305 
J = 27 P Mevls2T ies 2) 9.999993 %b 
J = 24 So = J-0103023)) 24 = 0.939999179 
J = 29 Ps Oeulo25277 2A = 9.99998989 
J = 30 Ps MEAT 71 FAs vA Ss 1.99992%759 
ys $1 SE Oe) 7 3 5 /m vA Ss 0.99998481 
Jos 32 PS 9.91772157 RAF 0.99998148 
J = 35 P= USO ARN1 313 eA s 0.99997749 
J = 34 P = OBUtR 30578) 9h = 0.9999727%3 
Ji s_ S5 Pos eT a> 2 S.ai oy = 0.9999670S 
J2= 36 Ps N21 4a 543 2A s 2.99996031 
J = 37 Pi GAO 7 a) QA s 0.999952 34 
Js 38 P= Je oo 7 OT QA = 9,9999429] 
j= os 9691975953 Q4 2 9.923993], AN 
Js an P= V.07A005910 a= 9.99991875 
er Ps 9.92033907 2A = 1,99990343 
a C= V.0705392 3 ee 0.99998SS1 
J = a3 Pos e120 72 1am PA = 0.99986459 
J = a4 Sy Oece tei 1s7 wis PL IIR) Ae 
Is YS P = ‘Peale Lasts vA = 9.99981189 
Jos 46 P= WeI2tI 725%) PA = 9.999779) 
Jos a7 as 0.072708 507 Se 4699973095 
Jos 48 » = 1.922 573n3 Qa = 029395969 
j = 49 P = Teese Syl e4 Ws os ¥.99454622 
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APPENDIX C 


COMPUTER PROGRAMS FOR CONVERTING 
FINITE ELEMENT (NIKE2D) PROGRAM OUTPUT TO BE 
STATISTICAL RELIABILITY COMPUTER PROGRAM INPUT 


Creare htkhhnhkh kh hr h dhe moe ett ek tr hetentkhe tenth etinhatkhhe eh 


AMrAMNMmAN 


PROGRAM: CONVERT NIKE20 OUTPUT TO BE STATISICAL INPUT 


FOR {PT. INTEGRATION ENTIRE SPECIMEN MODEL ONLY 


Creme trtkkae etek heheh heat hhhhnrehhhakheehhh 


OQ 


Creme kk RAR RARARARKERRAARARARRRARRAAAKRERERARRHEAKRaE ARKH AKKRAKh 


590 


100 


101 


103 


104 


105 


106 


107 


108 


OIMENSION XtL(594),xXwW(594), 
BYTE A(30) 
TYPE*,'ENTER FILE CONTAINING NIKE2D DATA QUTPUT' 
READ (§,500) A 
FORMATI30A1 ) 
OPEN(UNIT=S, FILE=A, STATUS='QLD’") 
OPENCUNIT=6, FILE='CC.OAT', STATUS='NEW' ) 
WRITE(6,50) 
NE=S9e 
NPSNE/2Z0+1 
DO 100 M=1,18 
¥L(M)=.1658334 
xW(M)=.095 
CONTINUE 
DO 101 M=19,36 
XL(M)2=.0995 
Xw(mM)=.095 
CONTINUE 
DO 102 M237,132 
XL(M)=.024875 
XW(M)=.095 
CONTINUE 
DO 193 M=333,150 
¥L(M)=,.0995 
Xwim)=.095 
CONTINUE 
DO 104 ™=1{51,168 
Xu (mM) = 21658334 
XW (M)=.095 
CONTINUE 
DO 105 M=169,200 
XL (M)=.1658334 
XW(M)=.024875 
CONTINUE 
DO 106 M=201,216 
xL(M)=,0995 
XW(M) 2.024875 
CONTINUE 
DO 107 M4=217,232 
XL (M)=.924875 
xw(M)=.928875 
CONTINUE 
DO 108 4=233,296 
XL(M)=,099S5 
XW(M) =.024875 
CONTINUE 
DO !09 M=297,316 
¥L OM) =.1658334 
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XN(M) =.024875 
109 CONTINUE 
DO 110 M=317,328 
XL (M) =. 1658334 
XW(M)=.024875 
1t0 CONTINUE 
DO 411{ M=329, 360 
XL (M)2.0995 
XW(M)=.024875 
i CONT INYE 
DOO 112 M=361,424 
X¥L(M)=.024875 
XW(M) 2.024875 
112 CONTINUE 
DO 113 M=42S, 456 
XL (4) 2.0995 
XW(M)=.024875 
113 CONTINUE 
DO 114 M=457,488 
XL (M)=. 1658334 
X4(M) 2.024875 
114 CONTINUE 
D0 115 4=4869,576 
X4(M) 5.024875 
xw(M12=.024875 
115 CONTINVE 
DG 116 4=577,592 
XL(M)=.024875 
XNW(M) 2.024875 
115 CONTINUE 
DO 19 oil »NP 
00 15 J=!,190 
M=lLOr( Jey) tu” 
TF (M.NE NE) THEN 
IF(S.EQ.10) THEN 
READ(5,45) SICY,SIGZ,SI3%,S(GY2Z 
ELSE - 
REAON(5,35) SIGY,SIGZ,SI135X,SIGYZ 
ENO IF 
ELSE 
READ(S5,30) SIGY,SIGZ,-S{GX,S16Y2 
60 TOeSs 
ENO [F 
WRITE(6,40) M,-XL(M),XWOM),SIGY,SIGZ,-SIGX,SIGYZ 
ES CONTINUE 
1G CONTINUE 
SS WRITE(6,40) M,XL(M) -XW(M),SIGY,SIGZ-SIG<,SIGYZ 
30 FORMAT(18X,4641.0) 
35 FORMAT(18X,4€611.0,////) 
4 FORMAT(3%,15,2F11.7,4€12.4) 
45 FORMAT(I8xX<,4uETT.O,///SISSSSISSSY 
50 FORMAT({X, ‘ELEMENT NO.',4X%,'L',10X,'W', 7X, 'SLGY',8x,'SIGZ' 
ep 8X. °STGX" 7X; sao 77) 
STOP 
END 


a 


Came t erate srt mi miniterrkehMethebanm eer mem ert hhtterreinkintenhheae 


CeCe ere rrr ree eee en ee ee ee ee a ee 2 2 2 2 2 2 2 2 2s 


Cc 
& PROGRAM: CONVERT NIKE20 ouTPuT TO BE STATISTICAL INPUT 
C 

C FOR 4pT. INTEGRATION HALF SPECIMEN MODEL ONLY 
Cc 


Crt meat r RAR RRkRRRARRREARARNERRERARH ERR ERERARKERKKRRERKHA 


Oo 


CRRA RRERRERARRE AAR ARERR RRARERRAARRRREERERRRARKRERARE 


DIMENSION XL( 1928) , XW( 4928) 
BYTE A(30) 
TyPEs, "ENTER FILE CONTALNING NIKE2D OQurPUT’ 
READ(S5,35090) A 
S00 FORMAT( 3041) 
OPEN(UNIT=S, FILESA, STATUSS'OLD') 
OPEN(UNTT=6, FILE='BB.04T', STATUS='NEW" ) 
WRITE (6,50) 
NE=482 
NEGINE RG 
NP=NE4/4Q + j 
00 100 M=1,{100 
X¥L(W)2.Q04975 
Xw(m)=.028590 
100 CONTINUE 
DO 404 M=101,220 
XL(M)=.024875 
XV(M)=.028500 
101 CONTINYE 
00 102 4=2411,540 
XL(M12.0124375 
xW0M)=.0285 000 
102 CONTINUE 
bo 493 m=541,660 
XL(4)5.024875 
x¥wim)=.028500 
103 CONTINUE 
00 104 M=661, 760 
XL() 2.04975 
xW(M)=.02850 
1904 CONTINUE 
DO 10S M=761,920 
XL(M)=.04975 
¥N(M)=.0124375 
10S “Gone 
DO 106 4=921,1412 
X£(M)=.024875 
XW(M)=.0124375 
106 CONFINUE 
00 107 M=3113,1368 
X¥L(M)=.0124375 
XW(M)=.0124375 
197 CONTINQE 
00 108 M=1369, 1560 
XL(M)=.024875 
xXW(m4)=.0124375 
108 COANTINYE 
DO 109 M=156),)1720 
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109 


110 


1 


20 
1S 


35 
39 


4Q 
a5 
20 


XL(M}=.0497S 
XW(M1=.0124375 
CONTINUE 
DG 110 mM=1721,1896 
XL(M?=.0124378 
XW(M) =.0124375 
CONTINUE 
00 111 M=1897,1928 
XL (M) =.0124375 
XW(M)=.0124375 
CONTINUE 
600 10 f=1,NP 
00 20 K=1,4 
MshOe(Lel) +4e( Jol) +K 
IF(M.NE.NEG) THEN 
IF(K.EQ.4) THEN 
IF (J.EQ.10) THEN 
READ(S,45) SIGY,SIGZ,SIGX,SIGYZ 
ELSE 
READ(S,35) SIGY,SIGZ,SIGX,SIGYZ 
END IF 
ELSE : 
READ(S5.30) SIGY,SIGZ,-SIG*,SIGYZ 
ENO [F 
ELSE 
REAO(S,30) SIGY,SIGZ,SIGX,.STGYZ 
GO TO 35S - 
END IF 
WRITE(6,40) MeXL CM) -XWOM),-SIGY-STGZ,SIGX,SLGYZ 
CONTINUE 
CONTINUE 
CONTINUE 
WRITE(6,40) M,XL CM), XW(M),STGY,SIGZ,SIGX,SISYZ 
FORMAT (48X,4E11.0) 
FORMAT(18x,4E611.0,7) 
FORMAT(3X,T5,2F11.7,-4€&12.4) 
FORMAT(18X%,4E14.0,////////7/) 
FORMAT(1%, "ELEMENT NO.',4%,'L',10K%,'W',7X,'SIGY', 8X, 'SIGZ' 
&,Qr, "SIGx' Pel or "SIrGYZ' +7 ) 
STOP 
END 


79 


Cette thea aReReke bibs et ReoRUAbDGRhCReRAhReCixPete Ve een tah se Pek oe 


- 


be 


Crt emake mekthaeraereh eae rekheonmehheaeeoeerheerteeheatenhtrekeeake 


PROGRAM: CONVERT NYJKE2D OUTPUT TO BE STATISTICAL INPUT 
FoR 4pT. INTEGRATION HALF SPECIMEN MODEL ONLY 


RkRkekartwarkkiknankkerenwrekaweeeterkavrkraekrreekekeaekRakanteeemnkre eee hee 


& 
C 
c 
E 
Carte keekaertketkakeaeeean ke ermkekhh Re herarherrtkeeaearhankhkeeeeea uk t 
G 
C 
e 


5090 


100 


191 


102 


103 


104 


105 


106 


107 


108 


DIMENSION XL( 482), xXW( 482) 
BYTE A(390) 
TYPE*, "FILE CONTAINING NIKE2D DATA QUTPUT' 
READ(S,500) A 
FORMAT (30A4) 
QOPENCUNIT=S, FILE=A, STATUS='OLO') 
OPEN(UNIT=6, FILE='B88B.DAT', STATUS='NEW' ) 
WRITE( 6,50) 
NE=482 
NP=NE/ 4041 
pd 1990 M=1,25 
XL(M)=.0995 
Kw(m)=.0570 
CONTINUE 
DO 101 M=26,55S 
XL (M) 2.04975 
XW(M)=.0S7 
CONTINUE 
DO 102 M=56,135 
XL(M)=.0249875 
XW(M)=.0570 
CONTINUE 
00 103 M=136,.16S 
XL(M)=.04975 
Xw(M) 2.0570 
CONTINUE 
DO 104 M=166,190 
XL(W)=.0995 
XW(M)=.0S70 
CONTINUE 
DO 1905 M=191,230 
XL(M) =.9995 
xw(mM) =.024875 
CONTINUE 
DO 106 M=231,.278 
XL(N)=.94975 
YW(M) =.024875 
CONTINUE 
06 107 M=279,342 
XL(M)=.024875 
XW(M)=.024875 
CONTINUE 
00 10A 4=343,390 
XL(M)=.04975 
XW(M)=.0246875 
CONTINUE 
DO 109 M=391,430 
XL(M)=.0995 
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XWiM) 2.024875 
199 CONTINUE 
DO 119 M2431,474 
Xu (M)=.024875 
XW(M)=.024875 
110 CONTINUE 
DO 111 M=475,482 
XL(M)=.024875 
XWEM) =.024875 
ated CONTINUE 
00 10 121,NP 
00 t§ J=1,19 
MsiOr(ITel)et 
IF (M.NE.NE) THEN 
IFC S.£EQ.10) THEN 
REAO(S,45) SIGY,SIGZ,SIGX,SIGYZ 
ELSE 
READ(5,35) SIGY,SIGZ,SIGX,SIGYZ 
ENO IF 
ELSE 
READ(S,39) SIGY,SIGZ,SIGX,SIGYZ 
GO TO $5 
ENO IF 
WRITE(6,40) Me¥L(M) -XW(M),SIGY,SIGZ,SIGX,SIGYZ 
iS CONTINUE 
10 CONTINUE 
SS WRITE(6,.40) MeXL(M) -XW(M),SIGY,SIGZ,SIGK,SIGYZ 
30 FORMAT (18xX,4E11.0) 
35 FORMAT(18XK,4E11.0-///7) 
1g FORMAT (3X ° 15,2eF 11 7 4E12. 4) 
YS FORMAT(14X,4E14 .0,7/4/7/77/7/77/7) 
S50 FORMAT(1%, ‘ELEMENT NO.',4AX,'L',JOX,'wW’,7x,'SIGY',3xX,'*SIGZ' 
#,8x,'SIGX%',7X,‘SIGYZ'.7/) 
STOP 
ENO 
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